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Abstract. This paper discusses some mathematical questions for spin-dependent Fermi gases at 
low temperatures modelled by a kinetic equation of Boltzmann type. The distribution functions 
have values in the space of positive definite hermitean 2x2 complex matrices. Global existence of 
c/3 ■ weak solutions to the initial value problem in n L°° for the Boltzmann equation is proved. 

ci ■ 

bJQ! 

1 Introduction. 

C ' 

^ ' The experimental study of spin polarized Fermi gases at low temperatures and their kinetic mod- 

elling is well established in physics, an early mathematical physics text in the area being [S]. The 
first experiments concerned very dilute solutions of "^He in superfluid '^He with - in comparison 
with classical Boltzmann gases - interesting new properites such as spin waves (see [NTLCL]). The 
experimentalists later turned to laser-trapped low temperature gases (see [JR]). The mathematical 
study of these models, however, is less advanced. To discuss this, we first recall some basics of the 
Pauli spin matrices 

1 \ f i \ f I 



^1 = 1 1 ; ' ^2 = ^ _. y , ^3 = ^ Q _^ 

Denoting by [cTi,cTj] the commutator a^aj — ajai, the Pauli matrices satisfy 

[o-i,cr2] = 2i(T3, [0-2,0-3] = 2icJi, [0-3,0-1] = 2io-2, and [o-j,o-j] = for i = l,2,3. (1.1) 

With 

a = (0-1,0-2,0-3) 
the Pauli spin vector, (1.1) is equivalent to o- x o- = 2io-. 

Let A^2(C) denote the space of 2 x 2 complex matrices with T-L2{C) the subspace of hermitean 
matrices. ?^2(C) is linearly isomorphic to using the decomposition p = A^I + As ■ a and 
identifying p G ^2(C) with {Ac,As) G M^. In the subspace of positive definite hermitean matrices, 
it holds 

log(p) = \{\og{Al - \AA') + \\og (^^^)^ • 

Zi Zi -^c 

The eigenvalues of logp are log(j4c — \ As\) and \og{Ac + \As\). For A,B ^ 7W2(C) the contracted 
product of A and B \s A : B = AijBij. The contracted product of two Pauli matrices o-j and aj 
gives o-j : aj = 26ij, where 6ij is Kronecker's delta. More generally for two vectors v = {vi,V2,V3) 
and w = {wi, W2,W3), it holds {v ■ a) : {w ■ a) = 2v ■ w. With / the identity matrix, {v ■ a){'w ■ a) = 
V • wl + i{v X w) ■ a, which implies [v ■ a,w • a] = 2i{y x w) • a. 
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A dilute spin polarized gas with spin ^, can be modelled kinetically by a distribution function 
matrix p with values in 7i2{C), which is the Wigner transform of the one-atom density operator for 
the system. In this paper p is assumed to evolve according to the kinetic equation 

Dp:=-p + p-sj^p = Q{p) (1.2) 
with in the Born approximation, the collision integral 

Q{p) = J dp2dp[dp'2B6{pi +P2- p'l - p'2)S{pI +pI- p'l - )({[pi, pi/]+TV(p2P2') 

-[Pl,Pl'] + Tr(P2P2')} - ^{[[^l,Pl'] + ,[P2,P2']+]+ - [[pi,Pl'],[p2,^2']+] + })- 

Here p = I — p, and [., .]+ denotes the anti-commutator. The domain of p{t,x,p) is positive time 
t, p G R^, and for simplicity periodic 3d position space x with period one. The kernel B is of 
hard force type, < B{\pi -p2\,d) = \pi - P2fb{e) with < P < 1, G L°°. The number 

density of particles is given by / := Tr(p(t, x, p)), and the magnetization of particles by the vector 
a{t,x,p) := Tr{ap(t,x,p)), which gives p = ^{fl + aa). The main result of the paper is a global 
existence theorem for the initial value problem of (1.2), 

Theorem 1.1 Suppose that (/o,(To) <= L°° D ^^([0, 1]-'^ x M^), and for p^ < f, j G N, that 
< r]j < fo < 2 — Tjj and o'q + r/| < min(/Q , (2 — /o)^) for some rjj > 0. Then the equation 
(1-2) with initial value (/o,^o); has a hounded integrable solution for t > with < / < 2 and 
a2<min(/2,(2-/)2) . 

Open problems of considerable physical interest are (cf [JM]) the relaxation times for spin-diffusion, 
the time asymptotic behaviour in general, and the influence of more involved transport terms in 
(1.2), such as the physicists' version of the problem, 

df ^ ^ sr^ \dh„ da dhp dai ^ 

't=xyz 

dt ^idpidri dri dpi dri dpi dpi (?r j J p m- 
Here Qn is the number density part and Qm the spin part of Q, the energy matrix is split into 

'v = ^ + j dp'{V{Q) - \v{\p~p'\)}nr,t), 

hp = -^{jB - J dp'V{p-p')a'ir,t), 

V is the inter-particle potential, B an external magnetic field, and 7 is the gyromagnetic ratio. 

Also linearized versions of (1.2) are discussed in [JM], but we are not aware of any related mathe- 
matical studies. However, related linear spinor Boltzmann equations introduced in spintronics for 
semiconductor hetero-structures, e.g. 
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■^P + v- VxP + E ■ \/yp = Q{p) + Qso{p) + Qsf{p), 



(1.4) 



have been analyzed mathematically. Here E is an electric field and Q is the collision operator for 
collisions without spin-reversal, in the linear BGK approximation 

/ a{v, v')iM{v)p{v') - M{v')p{v)dv' , 

with M denoting a normalized Maxwellian. The spin-orbit coupling generates an effective field 
making the spins precess. The corresponding spin-orbit interaction term Qso{p) is given by 
|[r2 • 0", p]. Finally Qsf{p) is a spin- flip collision operator, in relaxation time approximation given 

by 

Qsf{p) = , 

with Tgf > the spin relaxation time. Various mathematical properties - existence, uniqueness, 
asymptotic behaviour - have been studied (see [BH] and references therein). 



2 Preliminaries on the collision operator. 

The equations for / and a resulting from (1.2), are 

Df = QM,^), (2.1) 
Da = Qm{f,a), (2.2) 

where 

Qn(/i,^i) = dp2dp'idp'2B6{pi +P2-p'i -P2)Sipl +pl-p'i - p'i) 

(^([/l' - ^(/l/l' + ^1 • ^l')][f2' - ^(/2/2' + ^2 • ^2')] - [fl - ^{flfv + ^1 ■ ^l')][/2 - ^(/2/2' + ^2 ' ^2')]) 
-^([^1' - ^(/I'^l + /l^l')] • [^2' - ^(/2'0^2 + f2^2')] - [^1 - + fl^l')] ■ [^2 - ^(/2'0-2 + f2^2')])^ , 

Qm(/i, ^i) = ^ J dp2dp[dp'2B6{pi +P2- p'l - p'2)5{jpI +pI- p'l - p'i) 

(^([^1' - \Uv^l + h^v)] ■ [h - \U2f2' + ^2^2')] - [0^1 - ^(/I'^l + /l^l')][/2 - ^(/2/2' + ^2 • 0^2')]) 
-^([/l' - \{hh' + ^1 • ^V)][^2' - ^(/2'^2 + /2^2')] - [/l - \{hh' + ^1 • ^l')][^2 - ^(/2'^2 + /2^2')])) • 

The collision term does no change the number density (J Qndp = 0), the linear momentum 
density (J pQndp = 0), and the energy density (J p^Q^dp = 0), and the collisiion term does 

not change the magnetization density (J Qmdp = 0). 

Introduce for the direct sum of M and M^, the commutative ;-product given by (/i +a'i); (/2 + ^2) = 
/1/2 + cTi • CT2 + /ia-2 + f2^i- Obviously the ;-product is not in general associative, due to the inner 
product of a. Consider the sums Qn ± Qmi 
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Qn±Qm = ^ J dp2dp[dp2Bd{pi + P2 - p'l - p'2)^{p\ + pI - p'l - P2) 

([/2' - ^(/2/2' + ^2 • ^2')][(/l' ± ^1'); (1 - ^(/l ± ^1))] 
-[/2 - ^(/2/2' + 02 ■ a2')][{h ± ^1); (1 - \{fv ± aiO)] 

t[^2' - ^(/2'^2 + /20^2')]; [(/i' ± ^1'); (1 ^1))] 

±[^2 - \U2'&2 + /2C72'))]; [(/l ± ^1); (1 - \{fv ± ^1'))]) 

+\j dp2dp[dp'2Bdipi +P2- p'l - p'2)5ipl +pI- - p'i) 

([/2' - ^(/2/2' + ^2 • ^2')][(/l' ± ^1'); (1 - ^(/l ± ^1))] 

-[/2 - ^(/2/2' + 02 ■ a20][(/i ± ai); (1 - ± aiO)] 
±[^2' - ^(/2'^2 + 72^20]; [(/i' ± aiO; (1 - ^(/l ± ^1))] 

T[^2 - \U2'^2 + /2a2'))]; [(/l ± ^1); (1 - \{fv ± ay))]) 

= \j dp2dp[dp'2Bd{pi +P2- p'l - p'2)5{jpI +pI- p'l - p'i) 

([{h' ± ^1'); (1 - ^(/i ± ^1)]; [(/2' T ^2'); (1 - \U2 T ^2)] 

-[(/i ± ai); (1 - ± aiO]; [(/2 T ^2); (1 - ^(/2' T ^2')]) 
+\ J dp2dp'idp'2B5(jpi +P2- p'l - p'2)5{pl +pI- p'l - p'i) 

{[{fv ± av)- (1 - ^(/l ± ^1)]; [(/2' ± ^20; (1 - ^(/2 ± ^2)] 

-[(/i ± ai); (1 - \{fy ± ay)]- [{J2 ± ^2); (1 - \U2' ± ^2')]) • (2-3) 

In particular in the longitudinal case (all spins along the same axis) the equation reduces to the 
Uehling-Uhlenbeck equation for a spin-up, spin-down mixture, 

Dni± = ^J dpidp'idp'2S{pi +p2- p'l - p'2)S{pI +pI- Pi -p'h 

{^{rii^ri2±ni±n2± - ni±n2±rai±"2±) + ("i±'^2=F"'i±"2=F - '^i±?i2q=rai±n2=F)) , 
where n± = f ±a are the spin-up and spin-down densities. 
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For the product of (/i ^ ai) and (Q„ it Q^) it holds 

(/l - O^l); (<5n + Qm) + (/l + O^l); (<5n - Qm) = 

^ J dp2dp[dp'2Bd{pi +P2- p\ - p'2)5ipl +pI- p'^ - p'i) 

(/i - ^i); ([(/i' + ov)- (1 -\{h + ^i))]; [(/2' - ^2'); (i - \{h - ^2))]) 

-(/l - ^1); ([(/l + ^1); (1 - \{fv + [(/2 - ^2); (1 - \{f2' - &2'))]) 

+ih + ^1); ([(/r - aiO; (1 - ^(/i - ^1))]; [(/2' + ^2'); (1 - ^(/2 + ^2))]) 
-(/i + ai); ([(/i - ai); (1 - ^(/i, - ay))]; [(/a + ^2); (1 - ^(/2' + ^sO)]) + 
^ y dp2dp[dp'2Bd{pi +P2-p'i -P2)Sipl +pI-p'i-P2) 

ifi - ^1); {[{fv + aiO; (1 - ^(/i + ^1))]; [(/2' + ^2'); (1 - ^(/2 + ^2))]) 

-(/i - ai); ([(/i + ai); (1 - \ifr + a,,))]; [(/2 + ^2); (1 - ^(/2' + ^20)]) 
+ ^1); ([(/r - aiO; (1 - ^(/i - ^1))]; [(/2' - ^2'); (1 - ^(/2 - ^2))]) 
-(/i + ^1); ([(/i - ^1); (1 - \ih' - ^1'))]; [(/2 - ^2); (1 - ^(/2' - ^20)]). (2.4) 
The gain term part in the integrand of (2.4) is 

I ((/i - ^1); ([(/i' + ^1'); (1 - ^(/i + ^1))]; [(/2' - ^2'); (1 - ^(/2 - ^2))]) 

+(/i + ai); ([((/i. - ar); (1 - ^(/i - a,))]; [(/s^ + asO; (1 - ^(/2 + ^2))])) + 
J ((/i - ^1); ([(/i' + ^Ti^; : (1 - \{h + ai))]; [(/s' + ^2'); (1 - ^(/2 + ^^2))]) 
+(/i + ^1); ([(/i' - aiO; (1 - ^(/i - ^1))]; [(/2' - ^2'); (1 - lif2 - ^2))])) . 

(2.5) 

and the loss term part is 

\{-ih- ^1); ([(/i + ^1); (1 - \ih' + ^1'))]; [(/2 - ^2); (1 - \{f2' - ^2'))]) 

-(/l + ^1); ([(/l - ^1); (1 - lifl' - ^1'))]; [(/2 + ^2); (1 - ^(/2' + ^2'))])) + 

^ ( - (/i - ai); ([(/i + ai); (1 - ^{f,, + a,,))]; [(/2 + ^2); (1 - ^(/2' + a2^))]) 
-ifi + ai); ([(/i - ai); (1 - l(/r - aiO)]; [(/2 - ^2); (1 - ^(/2' - ^20)])). (2.6) 

In (2.5) the terms with one, three, or five cr-factors cancel out, since each a has opposite sign in the 
first and second line, and in the third and fourth line. The same holds for (2.6). 
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3 An existence result. 

We shall adapt an approach from [D] and [PLL] for Fermi-Dirac type Boltzmann equations, to study 
this initial value problem, and first treat the case with a truncation Bj in the domain of integration 
for Q, where Bj is the the restriction of B to the set Pi+p^ < i^- The initial value (/cO'o) £ 
for the equations (2.1-2) is assumed for each ball < f' to satisfy < r]j < fo < 2 — r]j, 
min((2 - /o)^, /q) > ctq • ctq + r/^ for some rfj > 0. Set 

F{t,x,p) = f{t,x,p) for 0</<2, =0 for / < 0, =2 for / > 2, E{t,x,p) = 

a{t,x,p) when min(F^, (2 - F)^) > a • a, else T.{t,x,p) = ^^^^^l^—^'^^ {t,x,p). 

ycr • a 

Consider the system 

Df = Qn{F,i:), Du = Qm{F,T.) (3.1) 
for t>Q under the truncation Bj and with the initial value (/o, ^o)(= (-^O; ^o))- 

Proposition 3.1 The system (3.1) with truncation Bj and initial value (/o,^o); has locally in time 
a unique solution (/, u) G L°° . 

Proof Set 

Tnf{t,x,p) = fo{x-tp,p)+ [ QniF,J:){s,x + {s-t)p,p)ds, (3.2) 

Jo 

Tma{t,x,p) = ao{x-tp,p)+ Qm{F,Y.){s,x + {s - t)p,p)ds. (3.3) 

^0 

Obviously (/, a) is a solution of the the initial value problem for (3.1), if and only if it is a fixed 
point of T := {Tn, T^). It is enough to prove that the operator T is contracting in L°° for < t < to, 
when < to is small enough. 
Now, 

r„(/l, ai) - r„(/2, C72) = ^j'dsj dp2dp[dp'2Bn5ipi + P2 - p[ - p'2)S{pl + pi - pf - p'i) 

(^2^[Fii' — -(FiiFii/ + Ell • 2ii/)][Fi2' — -{F12F12' + S12 • S12/)] 

— [-^11 — 2(-^ll-^ll' + ■ ^ll')][-^12 — 2iFl2Fi2' + S12 • Si2')]) 
-^(Pll' - ^(i^ll'Sii + FllSiiO] • [Si2' - ^(Fi2'Si2 + F12S12O] 

— Pii — 2^Fii'T,ii + FiiSii/)][Si2 — -(-Fi2'Si2 + -^2X12')])^ 

~\ J J '^P2#'ldP2-Bn^(pi +P2- p'l - P2)5(pl +pI- P? - P2) 
^y^^'i'V ^ -:^{F2\F2V + S21 • S21')][F22' — 2(-^22/22' + ^22 " S22')] 

— [F2I — 2 (-^21-^21' + S21 • S2l')][-^22 — 2 (-^22-^22' + ^22 • ^22')]) 
— 2(^21' — 2(-^21'S21 + -P21S21')] • [S22' — -^{F22i^22 + -^22^22')] 

— p21 — -(F21/S21 + F2iS21')][S22 — -(i*22'S22 + F22S22' )] ) 
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hence on [0, to] 

II Tn{Fi, El) — r„(F2, S2) \\l°°< toCniW Fi — F2 IIloo + II Si — S2 ||l°°) 

< ioC'ndl /l — /2 ||l°° + II ^1 — ^2 ||l°°)- 

An analogous estimate holds for T^n, and so T is contracting for to > and sufficiently small. □ 

It remains to prove that F = f and T, = a. This holds by continuity on a (short and j-dependcnt) 
time interval < t < tj, using the boundedness of T„ and T^- For the full interval [0, to] we first give 
the proof in the case when there is spin only in the as direction. With a = (0, 0, s), we can replace a 
with s in Qn and Qm, and the equation for a by the corresponding one for s. Notice in (2.4) that in 
the gain- and loss-parts separately, a spin-term will contain an odd number of ct's. But those terms 
have each different signs on the first and second, third and fourth lines of (2.5) and of (2.6), and so 
they cancel out. Terms with an even number of a's belong to M. In this spin-up spin-down case, the 
;-product is associative. Set ci = ISmaxp^ J dp2dp[dp2Bn5{pi +p2 — p'l — P2)S{p'i +P2 
It holds 

- Cliff - 4) < I dp2dp[dp'2B5{pi +P2- p[ - p'2)S{pI +pI- pf - p'i) 

l{-ih- si); ([(/i + ^1); (1 - ^(/i' + ^1'))]; [(/2 - ^2); (1 - ^(/2' - ^2'))]) 

-(/l + ^1); ([(/l - ^1); (1 - lifv - SV))]; [(/2 + S2); (1 - \{f2' + 52'))])) 

+\ ( - (/i - Si); ([(/i + Si); (1 - i(/r + sr))]; [(/2 + ^2); (1 - ^(/2' + 52'))]) 

-(/i + 5i); ([(/i - 5i); (1 - lifr - sy))]; [(/2 - ^2); (1 - ^(/2' - ^2'))])) 

< (/i - si); {Qn + Qm)(.f, s) + (/i + Si); {Qn - Qm){f, s) = D{fl - si). 

This implies 

< (/o - sl){x,p)e-'^'' < {f - sY{t,x,p), 

where (/, s) is the fixed point of (3.2-3) with initial value (/o,^o) • By continuity / it s > for 
< t < tj. Analogously, starting from the equation for (1 — ^{fi =b cri)) instead of /i it cri, we get 
by uniqueness the same solution to (3.2-3), together with the estimate 

< ((2 - /o)2 - sl){x,p)e-''^' < ((2 - /)2 - s^)*{t,x,p), (3.6) 

which by continuity holds for < t < tj together with 2 — / it s > for < t < tj. And so, 
uniformly in {x,p), \s\ < min(/, 2 — /) < 1 for < t < tj. By iteration, existence and uniqueness 
follow up to time to and then for t > 0. We conclude that p = \{fl + aa), solves the truncated 
initial value problem for (1.2) globally in time. 

The proof could also have been carried out with a instead of s, since the vector part of the in- 
tegrands is zero. Then by a slight extension, an analogous result follows in the case of general 
spin directions (general ao). Given {x,v), start from a spin base, where ai is in the aa-direction. If 
a[, a2, ^2 all arc in the 0-3-direction, then we are in the previous situation. Else with a = (o'l, 0^2, as), 
and a = (fTi,(j2) added, there are three additional terms with the new factors • CJ2, cr^ • cj2, and 
5-2 • 02. Since the spin-part of (/i - ai); {Qn + Qm){f, ^) + (/i + ^1); {Qn - Qm){f, ^) cancels, (3.4) 
still holds, now with is replaced by +1^1, and leading to 



(3.4) 
(3.5) 
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< ifi - ■ ao)(x,p)e-^i* < {f - ^ ■ a)*{t,x,p), 
< ((2 - fof - ao ■ ao)(x,p)e-'^i* < ((2 - ff - a ■ a)*{t,x,p). 

We conclude that p = \{fl + o'er) solves the truncated problem (3.1) in the general case, 

Proposition 3.2 The system (1-2) with truncation Bj and initial value {fo,^o), has a unique 
solution for t > 0. 

The initial value problem for (1.2) has a solution also in the limit p G M^. This follows using the 
weak compactness arguments for the Fermi-Dirac case of the Boltzmann equation in [PLL]. For 
this, we notice that the averaging lemma can be applied separately to each of the four component 
equations of (2.2). The entropy dissipation argument to control a restricted convergence of Q'j' 
in (134) of [PLL], can be replaced by the following direct control. Let K he a fixed compact set in 
Mp. For pi & K and j large, max(|p'^|, |p2l) > ^y^) when (B — Bj){pi — V2,uj) > 0. In the integral 
Q~^{pi,io) = J dp2 f d6\pi — P2\^b{6)f[f2{l — /i)(l — f2), make a change of variable from p2 to the 
one of Pi and p'2 giving that maximum, and estimate the other factors f and 1-f by one. With 
Ajipi,P2) = I ^tm'^ dO, this gives 

J dpi\Q+ - Q+\ = J dpidp2d9iB - Bj){pi -p2)/(/2(l - /i)(l - /2) 

<C J dp2f2 dpi{A - Aj){pi - P2)dpi ^0, j ^ 00. 

From here the proof in [PLL], now applied not just to the non-cancelling terms of /i/2(l — /i)(l~/2) 
and 71/2(1 — /()(1 — 72), but to all corresponding combinations of / and spin-components ak, k = 
1,2,3 holds, step by step. This proves Theorem 1.1. □ 
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